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ABSTRACT 


An  effective  computational  strategy  is  developed  for  generating  the  response  of  complex 
systems  using  (small  or  large)  perturbations  from  the  response  of  a  simple  structure  (or  a  simpler 
mathematical/discrete  model  of  the  original  structure).  Two  general  approaches  are  developed 
for  selecting  the  simpler  model  and  establishing  the  relations  between  the  original  and  simpler 
models.  The  two  approaches  are:  decomposition  or  partitioning  strategy,  and  hierarchical 
modeling  strategy.  Two  effective  partitioning  strategies  arc  used.  The  first  is  based  on  uncou¬ 
pling  of  load-carrying  mechanisms,  and  the  second  is  based  on  symmetry  transformations.  The 
hierarchical  modeling  used  is  a  predictor-corrector  iterational  process  based  on  using  a  simple 
mathematical  model  in  the  predictor  phase  and  correcting  the  response  using  a  more  accurate 
mathematical  model. 

RESEARCH  OBJECTIVES 

The  objective  of  the  present  study  is  to  develop  an  effective  computational  method  for 
generating  the  response  of  a  complex  system  using  large  perturbations  from  that  of  a  lower- 
order  model  associated  with  a  simpler  system  (or  a  simpler  mathematical/discrete  model  of  the 
original  system).  As  an  integral  part  of  the  proposed  strategy  an  attempt  will  be  made  to  unify 
and  realize  the  full  potential  of  a  number  of  multilevel  computational  strategies,  some  of  which 
were  developed  by  the  principal  investigator  and  his  colleagues.  The  multilevel  strategies 
include  reduction  methods,  hybrid  modeling! analysis  techniques,  and  partitioning  methods. 
Reduction  methods  are  techniques  for  substantially  reducing  the  number  ol  degrees  of  freedom  of 
the  initial  discretization,  and  have  been  successfully  applied  to  a  number  of  vib  ution  and  non¬ 
linear  structural  and  thermal  problems.  Hybrid  modeling/analysis  techniques  can  achieve 
significant  reductions  in  the  analysis  time  by  incorporating  the  known  physical  behavior  into  the 
computational  model  of  the  system  and  by  using  different  analysis  methods  and/or  models  in 
predicting  the  different  response  characteristics  of  the  engineering  systems.  Partitioning  meth¬ 
ods  arc  bused  on  bieaking  the  large  (and/or  complex)  problem  into  a  number  of  smaller  (and/or 
simpler)  subproblems.  The  solution  of  the  original  problem  is  generated  using  information 
provided  by  the  individual  subproblem:: 


The  proposed  strategy  is  believed  to  combine  the  following  three  major  characteristics: 

1 )  gives  physical  insight  about  the  response 

2)  helps  in  assessing  the  adequacy  of  the  computational  model;  and 

3)  is  highly  efficient. 

The  strategy  will  first  be  applied  to:  a)  the  nonlinear  postbuckling  problem  of  composite 
structures ;  b)  reanalysis  of  large  structures  in  the  presence  of  geometric  nonlineariiies ;  then  c) 
coupled  field  problems.  The  postbudding  response  of  die  highly  anisotropic  composite  structure 
is  generated  using  large  perturbations  from  the  response  of  a  simpler  structure.  The  three  key 
elements  of  the  strategy  to  be  exploited  in  the  first  two  applications  are:  1)  mixed  (or  primitive 
variable)  formulation,  with  the  fundamental  unknowns  consisting  of  both  stress  and  displacement 
parameters;  2)  operator  splitting,  or  additive  decomposition  of  the  different  arrays  in  the  equa¬ 
tions  of  the  given  structure  to  the  corresponding  arrays  of  the  simpler,  or  previously-analyzed, 
structure  plus  correction  terms;  and  c)  application  of  a  reduction  method  and/or  a  stable  iterative 
method  for  the  efficient  generation  of  the  equations  of  the  given  structure. 

RESEARCH  ACCOMPLISHMENTS 

During  the  period  March  1,  1990  to  August  31,  1990,  two  tasks  have  been  performed.  The 
first  is  the  development  of  an  improved  partitioning  strategy  for  large-scale  structural  problems. 
The  second  is  the  development  of  a  predictor-corrector  approach  for  generating  the  steady-state 
thermal  response  of  multilayered  composite  plates  and  shells.  The  two  tasks  are  described 
subsequently. 

Improved  Partitioning  Strategy  for  Large-Scale  Problems 

The  governing  equations  for  the  discrete  model  of  the  original  structure  can  be  written  in 
the  following  compact  form: 

|K||Z)  =  {Q)  (1) 

where  { Z )  is  the  vector  of  stress  parameters  and  generalized  displacements;  |K|  is  the  global 
structure  matrix  which  includes  the  flexibility  and  strain-displacement  matrices;  and  (Q)  is  the 
global  right-hand-side  vector. 

In  the  decomposition  strategy  the  vector  of  fundamental  unknowns  { Z I  is  partitioned  into 


smaller  subvectors.  The  governing  discrete  equations,  Eqs.  1,  are  partitioned  accordingly.  The 
simpler  model  is  associated  wih  the  uncoupled  equations  in  the  partitioned  variables.  For  the 
case  of  two  partitions,  the  process  can  be  described  by  embedding  Eqs.  I  in  a  single-parameter 
family  of  equations  as  follows: 


(2) 


where  [Z| } ,  {Z2}  and  {Q[ } ,  {Q2}  are  the  partitions  of  the  original  vectors  { Z }  and  (QJ;  X  is  a 
tracing  parameter  which  identifies  all  the  correction  terms  needed  in  going  from  the  simpler 


model  to  the  discrete  model  of  the  original  structure;  and  a  dot  (.)  refers  to  a  zero  submatrix.  The 


simpler  model  corresponds  to  X=()  (uncoupled  equations  in  { Z 5 }  and  {Z2}),  and  the  discrete 
model  of  the  original  structure  corresponds  to  X=\  (fully  coupled  equations).  The  solution 
corresponding  to  A.=  l  is  generated  from  the  corresponding  solution  at  A.=()  using  an  iterative 
process  such  as  the  Preconditioned  Conjugate  Gradient  (PCG).  Note  that  the  correction  vectors 
of  the  iterative  process,  provide  a  direct  measure  of  the  sensitivity  of  the  response  quantities  to 
the  coupling  terms  (viz.,  the  terms  associated  with  the  tracing  parameter  X  in  Eq.  2). 

The  vectors  {Z||  and  {Z2}  are  chosen  to  be  the  symmetric  and  antisymmetric  compo¬ 
nents  of  the  response  vector  (each  is  approximately  half  the  size  of  the  original  vector,  { Z } ).  The 
simpler  model  (A.=0)  corresponds  to  a  symmetrized  structure  in  which  the  symmetric  and  an¬ 
tisymmetric  components  of  the  response  vector  are  uncoupled.  This  approach  can  be  thought  of 
as  a  physical  domain  decomposition.  If  the  PCG  technique  is  used  in  generating  the  solution  at 
A  =  l,  and  the  preconditioning  matrix  is  selected  to  be  the  left-hand  side  matrix  corresponding  to 
A=(),  then  each  of  the  correction  vectors  is  either  symmetric  or  antisymmetric. 

The  convergence  ol  the  PCG  technique  can  be  expedited  by  replacing  Eqs.  2  by  the 
following  equivalent  form  of  two  uncoupled  equations  in  { Z,  \  and  { Z2 }: 


(3) 


Note  that  for  X=l,  each  diagonal  block  of  the  total  left-hand-side  matrix  is  in  the  form  of  Schur 


complement  which  is  not  formed  explicitly.  Rather,  the  preconditioning  matrix  is  selected  to  lx* 

the  first  matrix  on  the  left-hand  side  (corresponding  to  X=0)  and  the  PCG  technique  is  used  in 

generating  the  solution  at  A.=  l .  The  results  of  this  research  are  reported  in  Ref.  1 . 

Predictor-Corrector  Approach  for  Generating  the  Steady-State  Thermal  Response  of  Multi¬ 
layered  Plates  and  Cylinders 

A  predictor-corrector  procedure  has  been  developed  for  the  accurate  determination  of  the 
temperature  and  heat  flux  distributions  in  thick  multilayered  composite  plates  and  shells.  The 
procedure  is  based  on  using  a  linear  through-the-thickness  temperature  distribution  in  the  predic¬ 
tor  phase.  The  functional  dependence  of  temperature  on  the  thickness  coordinate  is  then  calcu¬ 
lated  a  posteriori  and  used  in  the  corrector  phase. 

Extensive  numerical  results  have  been  conducted  for  linear  steady-state  heat  conduction 
problems,  showing  the  effects  of  variation  in  the  geometric  and  lamination  parameters  <  i  the 
accuracy  of  the  thermal  response  predictions  of  the  predictor-corrector  approach.  Both  antisym¬ 
metrically  laminated  anisotropic  plates,  and  multilayered  orthotropic  cylinders  are  considered. 
The  solutions  are  assumed  to  be  periodic  in  the  surface  cordinates.  For  each  problem  the  stand¬ 
ard  of  comparison  is  taken  to  be  the  analytic  three-dimensional  solution  based  on  treating  each 
layer  as  a  homogeneous  anisotropic  medium.  The  potential  of  the  predictor-corrector  approach 
for  predicting  the  thermal  response  of  multilayered  plates  and  shells  with  complicated  geometry 
is  investigated.  The  results  of  this  study  are  reported  in  Ref.  2. 

PUBLICATIONS 

1.  Noor,  A.  K.  and  Peters,  J.  M.,  "Strategies  for  Large-Scale  Structural  Problems  on  High- 
Performance  Computers,"  Communications  in  Applied  Numerical  Methods  (to  appear). 

2.  Noor,  A.  K.,  Steady-State  Heat  Conduction  in  Multilayered  Composite  Plates  and  Shells," 
Computers  and  Structures  (to  appear). 
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SUMMARY 

Novel  computational  strategies  are  presented  for  the  analysis  of  large  and  complex  struc¬ 
tures.  The  strategies  are  based  on  generating  the  response  of  the  complex  structure  using  large 
perturbations  from  that  of  a  simpler  model,  associated  with  a  simpler  structure  (or  a  simpler 
mathematical/discrete  model  of  the  original  structure).  Numerical  examples  are  presented  to 
demonstrate  the  effectiveness  of  the  strategies  developed. 


STEADY-STATE  HEAT  CONDUCTION  IN  MULTILAYERED 
COMPOSITE  PLATES  AND  SHELLS 
Ahmed  K.  Noor  and  W.  Scott  Burton 
Center  for  Computational  Structures  Technology 
NASA  Langley  Research  Center 
Hampton,  VA  23665 

ABSTRACT 

A  study  is  made  of  a  predictor-corrector  procedure  for  the  accurate  determination  ol  the 
temperature  and  heat  flux  distributions  in  thick  multilayered  composite  plates  and  shells.  A 
linear  through-the-thickness  temperature  distribution  is  used  in  the  predictor  phase.  The  lunc- 
tional  dependence  of  temperature  on  the  thickness  coordinate  is  then  calculated  a  posteriori  and 
used  in  the  corrector  phase. 

Extensive  numerical  results  are  presented,  for  linear  steady-state  heat  conduction 
problems,  showing  the  effects  of  variation  in  the  geometric  and  lamination  parameters  on  the 
accuracy  of  the  themial  response  predictions  of  the  predictor-corrector  approach.  Both  antisym¬ 
metrically  laminated  anisotropic  plates,  and  multilayered  orthotropic  cylinders  are  considered. 
The  solutions  are  assumed  to  be  periodic  in  the  surface  coordinates.  For  each  problem  the 
standard  of  comparison  is  taken  to  be  the  analytic  three-dimensional  solution  based  on  treating 
each  layer  as  a  homogeneous  anisotropic  medium.  The  potential  of  the  predictor-corrector 
approach  for  predicting  the  thermal  response  of  multilayered  plates  and  shells  with  complicated 
geometry  is  discussed. 
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STRATEGIES  EOR  LARGE-SCALE  STRUCTURAL  PROBLEMS 
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Ahmed  K.  Noor  and  Jeanne  M.  Peters 
Center  for  Computational  Structures  Technology 
University  of  Virginia 
NASA  Langley  Research  Center 
Hampton,  VA  23665 

SUMMARY 

Novel  computational  strategies  are  presented  for  the  analysis  of  large  and  complex  struc¬ 
tures.  The  strategies  are  based  on  generating  the  response  of  the  complex  structure  using  lar^c 
perturbations  from  that  of  a  simpler  model,  associated  with  a  simpler  structure  (or  a  simpler 
mathematical/discrete  model  of  the  original  structure).  Numerical  examples  are  presented  to 
demonstrate  the  effectiveness  of  the  strategies  developed. 

1.  INTRODUCTION 

Dynamic  and  nonlinear  response  calculations  for  future  aerospace  systems  are  expected  to 
require  processing  rates  far  in  excess  of  computers  built  around  a  single  processing  unit.  This  is 
because  of  the  complexity  of  these  systems  and  the  high  degree  of  sophistication  of  the  comput¬ 
ational  models  required  for  simulating  their  response  (see,  for  example.  Refs.  I  and  2).  The 
introduction  of  novel  forms  of  machine  architecture  has  brightened  the  prospects  for  meeting 
future  large-scale  computational  needs.  Most  of  the  new  machines  achieve  high  performance 
through  vectorization  and/or  parallelism.  These  include  top-of-the-range  supercomputers  such  as 
CRAY-3,  CRAY  Y-MP,  as  well  as  computing  systems  based  on  readily  available  and  inexpen¬ 
sive  basic  units  such  as  the  hypercubes,  the  Connection  Machine,  and  the  transputer  networks. 
The  characteristics  of  several  new  machines  are  summarized  in  Refs.  3,  4,  5  and  6.  Much  work 
has  been  devoted  to  the  development  of  efficient  vector  and  parallel  numerical  algorithms  lor 

r 

performing  matrix  operations,  solution  of  algebraic  equations,  and  extraction  of  eigenvalues  (see. 
for  example.  Refs.  7,  8  and  9).  Also,  a  number  of  special  strategies  have  been  proposed  for 
increasing  the  degree  of  parallelism  and/or  vectorization  in  finite  element  computations.  Most  of 
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the  special  strategies  can  be  thought  of  as  applications  of  multilevel  computational  processes,  and 
include  reduction  methods,  hybrid  modeling/analysis  techniques,  and  partitioning  methods. 
Reduction  methods  arc  techniques  for  reducing  the  number  of  degrees  of  freedom  of  the  initial 
discretization  and  have  been  successfully  applied  to  a  number  of  vibration  and  nonlinear  prob¬ 
lems  of  structures  (see,  for  example,  Refs.  10-13).  Hybrid  modelingl analysis  techniques  can 
achieve  significant  reductions  in  the  analysis  time  by  incorporating  the  known  physical  behavior 
into  the  computational  model  of  the  structure  and  by  using  different  analysis  methods  and/or 
models  in  predicting  the  different  response  characteristics  of  the  structure  (see.  for  example, 
Refs.  14  and  15).  Partitioning  methods  are  based  on  the  intuitively  obvious  and  well-established 
practice  of  breaking  the  large  (and/or  complex)  problem  into  a  number  of  smaller  (and/or 
simpler)  subproblcms.  The  solution  of  the  original  problem  is  generated  using  information 
provided  by  the  individual  subproblems  (see,  for  example,  Ref.  16).  The  two  key  advantages  of 
partitioning  techniques  are  computational 
efficiency  and  modular  implementation. 

The  present  study  is  an  attempt  to  unify  and  realize  the  full  potential  of  a  number  of 
multilevel  computational  strategies  for  solution  of  large-scale  structural  problems.  Specifically, 
the  objective  of  the  present  paper  is  to  report  the  progress  made  in  developing  computational 
strategies  which,  in  addition  to  being  efficient  on  high-performance  computers,  combine  the 
following  two  major  characteristics: 

1)  give  physical  insight  about  the  response;  and 

2)  help  in  assessing  the  adequacy  of  the  computational  model. 

2.  BASIC  IDEA  AND  KEY  ELEMENTS  OE  THE  PROPOSED  STRATEGIES 

The  strategies  developed  are  based  on  generating  the  response  of  a  large  and  complex 
structure  using  larye  perturbations  from  that  of  a  simpler  model  associated  with  a  simpler 
structure  (or  a  simpler  mathematical/discrete  model  of  the  original  structure).  The  key  elements 
of  the  proposed  strategies  tire:  1)  mixed  (primitive  variable)  formulation  with  the  fundamental 
unknowns  consisting  of  stress  parameters,  generalized  displacements  and,  for  dynamic  problems, 
velocity  components;  2)  operator  splitting  or  restructuring  of  the  discrete  equations  of  the  origi- 
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mil  (complex)  structure  to  delineate  the  contributions  of  the  simpler  model  ;'.nd  the  correction 
terms;  and  3)  stable  iterative  process  for  the  efficient  generation  of  the  response  of  the  complex 
structure  starting  from  that  of  the  simpler  structure.  Although  in  some  applications  the  strategies 
can  be  used  in  conjunction  with  the  single-field  displacement  formulation,  there  tire  several 
applications  in  which  the  use  of  the  mixed  formulation  is  essential. 

2  i  Cjovembig  Equations 

The  governing  equations  for  the  discrete  model  of  the  original  structure  can  be  written  in 
the  following  compact  form: 

!K1(Z)  =  (0)  (I) 

where  (Z)  is  the  vector  of  stress  parameters  and  generalized  displacements;  |K|  is  the  global 
structure  matrix  which  includes  the  flexibility  and  strain-displacement  matrices;  and  (0)  is  the 
global  right-hand-side  vector. 

2.2  Relations  Between  the  Original  and  Simpler  Models 

The  crux  of  the  proposed  strategies  is  the  proper  selection  of  the  simpler  structure  and/or 
the  simpler  model,  for  a  given  complex  structure,  a  multitude  of  choices  of  the  simpler  structure 
can  be  made.  The  best  choice  is  that  which  exploits  the  characteristics  of  the  high-performance 
computers  and  satisfies,  to  the  greatest  extent  possible,  the  two  major  characteristics  lister!  in  the 
preceding  section,  for  each  choice,  the  relations  between  the  simpler  model  and  the  original 
structure  can  be  established,  and  an  efficient  computational  procedure  can  be  developed  for 
generating  the  response  of  the  original  structure  using  large  perturbations  from  that  of  the 
simpler  model.  Two  general  approaches  arc  described  subsequently  for  selecting  'he  simpler 
model  and  establishing  the  relations  between  the  original  and  simpler  models.  The  two  ap¬ 
proaches  are:  decomposition  or  partitioning  strategy,  and  hierarchical  modeling  strategy.  For 
simplicity,  the  two  approaches  are  described  with  reference  to  their  application  to  linear  stress 
analysis  problems.  I  lowever,  they  are  equally  applicable  to  nonlinear  and  dynamic  problems. 

l' 

3.  DECOMPOSITION  OR  PARTITIONING  STRATEGY 

In  the  decomposition  strategy  the  vector  of  fundamental  unknowns  (Z|  is  partitioned  into 
smaller  subvectors.  The  governing  discrete  equations,  Eqs.  I,  are  partitioned  accordingly.  The 
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simpler  model  is  associated  will  the  uncoupled  equations  in  the  partitioned  variables.  For  the 
case  of  two  partitions,  the  process  can  lie  described  by  embedding  Fqs.  1  in  a  single-parameter 
family  of  equations  as  follows: 


(2) 


where  {Z\} ,  {Z2}  and  {Qi  }  ,  {Q2}  are  die  partitions  of  the  original  vectors  (Z)  and  {01;  X  is  a 
tracing  parameter  which  identifies  all  the  correction  terms  needed  in  going  from  the  simpler 

model  to  the  discrete  model  of  the  original  structure;  and  a  dot  (.)  refers  to  a  zero  submatrix.  The 
simpler  model  corresponds  to  X=0  (uncoupled  equations  in  {Z|}  and  {Z2}),  and  the  discrete 
model  of  the  original  structure  corresponds  to  X=1  (fully  coupled  equations).  The  solution 
corresponding  to  X=  1  is  generated  from  the  correspond!'  •  solution  at  X=0  using  an  iterative 
process  such  as  the  Preconditioned  Conjugate  Gradient  (PCG).  Note  that  the  correction  vectors 
of  the  iterative  process,  provide  a  direct  measure  of  the  sensitivity  of  the  response  quantities  to 
the  coupling  terms  (viz.,  the  terms  associated  with  the  tracing  parameter  X  in  Eq.  2).  Among  the 
effective  partitioning  strategies  developed  in  the  present  study  are  the  following  two: 

a)  Uncoupling  of  Load-Carrying  Mechanisms.  For  example,  if  the  original  structure  is 
associated  with  a  two-dimensional  shell  structure,  the  simpler  structure  can  be  associated  with 
the  corresponding  plate  structure  in  which  the  membrane  and  bending  load-carrying  mechanisms 
are  uncoupled. 

b)  Symmetry  Transformations.  This  is  accomplished  by  selecting  {Z|}  and  { Zi }  to 
represent  the  symmetric  and  antisymmetric  components  of  the  response  vector  (each  is  ap¬ 
proximately  half  the  size  of  the  original  vector,  (Zj).  The  simpler  model  (X=0)  corresponds  to  a 
symmetrized  structure  in  which  the  symmetric  and  antisymmetric  components  of  the  response 
vector  are  uncoupled.  This  approach  can  be  thought  of  as  a  physical  domain  decomposition.  If 
the  PCG  technique  is  used  in  generating  the  solution  at  X=J,  and  the  preconditioning  matrix  is 
selected  to  be  the  left-hand  side  matrix  corresponding  to  X=(),  then  each  of  the  correction  vectors 
is  cither  symmetric  or  antisymmetric.  The  details  of  the  symmetry  transformations  are  described 
in  Refs.  17,  IK  and  20. 
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4.  COMMENTS  ON  THE  DECOMPOSITION  STRATEGY 

The  following  comments  regarding  the  decomposition  strategy  are  in  order: 

1.  The  convergence  of  the  PCG  technique  can  be  expedited  by  replacing  Eqs.  2  by  the 
following  equivalent  form  of  two  uncoupled  equations  in  |Z| }  and  { Zo  } : 

•  ]]f'LH-4K,2K“02l  «j, 

K21  Ki!  K12jJlZ2j  1°2J  |K21KuQ,J 

Eqs.  3  are  obtained  from  Eqs.  2  by  eliminating  {Z2}  from  the  first  equation  and  { Z|  }  from  the 
second  equation.  Note  that  for  X=l,  each  diagonal  block  of  the  total  left-hand-side  matrix  is  in 
the  form  of  Schur  complement  which  is  not  formed  explicitly.  Rather,  the  preconditioning 
matrix  is  selected  to  be  the  first  matrix  on  the  left-hand  side  (corresponding  to  X=0)  and  the  PCG 
technique  is  used  in  generating  the  solution  at  A.=  l.  The  technique  is  outlined  in  Appendix  A. 

2.  Equations  3  can  be  thought  of  as  the  governing  finite  element  equations  of  an  equiva¬ 
lent  structure  whose  response  under  the  modified  loads  (given  by  the  total  right-hand-side 
vector ).  is  identical  to  the  response  of  the  actual  structure  when  subjected  to  the  given  loads. 
The  equivalent  structure  has  uncoupled  load-carrying  mechanisms  (or  uncoupled  symmetric/ 
antisymmetric  response  vectors). 

3.  The  decomposition  process  can  be  repeated  to  effect  further  reduction  in  the  size  of  the 
partitions  of  the  original  vectors  and  matrices.  However,  us  the  number  of  partitions  increases, 
the  coupling  between  them  increases,  and  therefore,  the  number  of  iterations  in  the  PCG  tech¬ 
nique  increases.  Consequently,  the  proposed  strategy  may  not  be  effective  on  multiprocessor 
computers  with  fine  granularity  and  small  local  memories  (e.g.,  the  hypercubes  and  the  Connec¬ 
tion  Machine  of  Thinking  Machines,  Inc.). 

5.  HIERARCHICAL  MODELING  STRATEGY 

/* 

The  simpler  model  is  selected  as  one  with  considerably  fewer  degrees  of  freedom  than 
those  of  the  discrete  model  of  the  original  structure.  The  governing  equations  of  the  simpler 
model  can  be  written  in  an  analogous  form  to  those  of  the  original  structure,  as  follows: 
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|k|(/.)  =  (q)  (4) 

Tlie  response  vectors  of  the  original  and  simpler  models  are  related  through  the  interpolation 
operator  (F|  as  follows: 

IZ)  =  |r|(z)  (5) 

where  |T|  is  a  rectangular  transformation  matrix  ( interpolation  operator). 

Based  on  IZqs.  4,  the  relations  between  fkj,  |KJ  and  (q),  |Q|  are  given  by: 

ikj  =  in'iKjin  (6) 

lq)=ir|'{Q)  (7) 


where  superscript  t  denotes  transposition,  and  [T]1  is  referred  to  as  the  restriction  operator. 

Among  the  possible  choices  of  the  simpler  model  in  this  category  are  the  following  two: 

a)  Discrete  model,  based  on  the  same  mathematical  model  and  the  same  discretization 
procedure  used  for  the  original  structure,  but  with  considerably  fewer  degrees  of  freedom.  The 
resulting  numerical  process  is  similar  to  the  classical  multigrid  technique. 

b)  Discrete  model,  based  on  a  mathematical  model  of  lower  dimensionality  than  that  used 
for  the  original  structure.  For  example,  if  the  original  structure  is  modeled  using  a  two- 
dimensional  plate  or  shell  theory,  the  simpler  model  can  be  associated  with  a  one-dimensional 
thin-walled  beam  theory  (with  the  effects  of  flexural-torsional  coupling  included).  The  interpola¬ 
tion  operator,  IT],  then  reflects  the  basic  assumptions  made  in  the  dimensionality  reduction  (i.e„ 
in  going  from  a  two-dimensional  shell/plate  theory  to  a  one-dimensional  thin-walled  beam 
theory;  namely,  the  projection  of  each  cross-section  on  a  plane  normul  to  the  initial  cent roiclal 
axes  does  not  distort  during  deformation).  This  approach  can,  therefore,  be  thought  of  as  a 
physical  multigrid  method. 


6.  NUMERICAL  STUDIES 

To  assess  the  effectiveness  of  the  foregoing  strategics*  a  number  of  linear  and  nonlinear 
stress  analysis,  free  vibration,  buckling  and  dynamic  problems  of  complex  structures  have  been 
solved  by  these  strategies.  For  each  problem  the  solutions  obtained  by  the  foregoing  strategies 
were  compared  with  those  obtained  by  the  direct  analysis  of  the  original  (complex)  structure. 
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Some  of  these  applications  are  reported  in  Refs.  17-20.  Two  applications  of  the  decomposition 
strategies  are  presented  herein:  1)  linear  stress  analysis  of  cantilevered  composite  shallow  shell 
with  unsymmetrical  lamination  in  the  thickness  direction;  2)  nonlinear  dynamic  problem  of 
laminated  anisotropic  panel  with  an  off-center  circular  cutout.  The  two  applications  are  dis¬ 
cussed  subsequently. 

Both  structures  were  modeled  by  using  mixed  finite  element  models  with  the  stress  resul¬ 
tants  allowed  to  be  discontinuous  at  interelement  boundaries.  A  12x12  grid  was  used  for  the  first 
structure  and  a  192-element-grid  was  used  for  the  second  structure.  Biquadratic  shape  functions 
were  used  in  approximating  each  of  the  generalized  displacements,  and  bilinear  shape  functions 
were  used  for  approximating  each  of  the  stress  resultants.  The  characteristics  of  the  finite 
element  model  are  given  in  Ref.  21. 

6.1  Linear  Stress  Analysis  of  Cantilevered  Composite  Shallow  Shell 

As  an  application  of  the  strategies  based  on  uncoupling  of  the  load-carrying  mechanisms 
and  symmetry  transformations,  consider  the  graphite-epoxy  composite  shallow  shell  with 
trapezoidal  planform  shown  in  Fig.  1.  The  shell  is  subjected  to  uniform  normal  loading  and  has 
nonzero  curvatures  k,  and  k2  and  an  unsymmetric  lamination  in  the  thickness  direction.  The 

bending  and  membrane  load-carrying  mechanisms  of  the  structure  are,  therefore,  coupled.  Also, 
because  of  the  unsymmetry  of  the  structure,  the  response  does  not  exhibit  any  symmetry.  The 
two  decomposition  strategies  outlined  in  section  3  were  applied.  The  number  of  stress  and 
displacement  degrees  of  freedom  in  each  partition  is  given  in  Table  1. 

a)  Uncoupling  of  Load-Carrying  Mechanisms.  The  simpler  structure  is  selected  to  be  the 
composite  plate  with  zero  curvatures  and  with  the  bending-extensional  coupling  neglected 
(corresponding  to  X=0  in  Eqs.  2  and  3). 

Figure  2  shows  the  normalized  contour  plots  for  the  generalized  displacements  in  the 
original  and  simpler  structures.  The  strong  coupling  between  the  in-plane  and  bending  degrees 
of  freedom  is  evident  from  this  figure.  An  indication  of  the  accuracy  and  convergence  of  the 
solutions  obtained  by  using  the  PCG  technique,  in  conjunction  with  both  Eqs.  2  and  3,  is  given  in 
Fig.  3.  The  standard  of  comparison  is  taken  to  be  the  direct  solution  of  the  original  structure.  As 
can  be  seen  from  Fig.  3,  the  normal  displacement  and  total  strain  energy  obtained  by  using  the 
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PCG  technique,  in  conjunction  with  both  Eqs.  2  and  3,  converge.  However,  the  convergence  is 
faster  for  Eqs.  3. 

b)  Symmetry  Transformations.  The  simpler  structure  is  selected  to  be  a  structure  whose 
response  vector  exhibits  mirror  symmetry  with  respect  to  line  cc  (see  Fig.  4).  Solutions  were 
obtained  using  the  PCG  technique,  in  conjunction  with  Eqs.  2  and  3.  Also,  for  the  sake  of 
comparison  the  PCG  technique  was  used  in  conjunction  with  the  displacement  finite  element 
model  which  is  equivalent  to  the  mixed  model  used  herein  (see  Ref.  21).  The  vecluis  i Z, )  and 

(Zn)  in  Eqs.  2  and  3  arc  associated  with  displacement  degrees  of  freedom  only.  Figure  4  shows 
the  normalized  contour  plots  for  the  displacements  in  the  original  and  simpler  structures.  An 
indication  of  the  accuracy  and  convergence  of  the  solutions  obtained  by  using  the  different  PCG 
solutions  is  given  in  Fig.  5.  As  can  be  seen  from  Fig.  5,  the  displacement  formulation  results  in 
much  slower  convergence  than  the  mixed  formulation  and  for  both  the  mixed  and  displacement 
formulations,  the  use  of  Eqs.  3  results  in  faster  convergence  than  that  of  Eqs.  2. 

6.2  Nonlinear  Dynamic  Analysis  of  an  Anisotropic  Panel  with  an  Off-Center  Circular  Cutout 

As  an  application  of  the  strategy  based  on  symmetry  transformations  to  nonlinear  dynamic 
problems,  consider  the  laminated  anisotropic  cylindrical  panel  with  an  off-center  circular  cutout 
shown  in  Fig.  6.  The  loading  is  assumed  to  be  uniformly  distributed  and  normal  to  the  panel 
surface  and  has  a  step  variation  in  time.  The  panel  is  made  of  graphite-epoxy  material  (see  Fig. 
6).  Mixed  finite  element  models  were  used  for  the  spatial  discretization,  and  implicit  three-step 
method  was  used  for  the  temporal  integration. 

A  two-level  iterative  process  (with  two  nested  iteration  loops)  was  used  to  generate  the 
symmetric/antisymmetric  components  of  the  response  vectors  at  each  time  step.  The  top  level 
iteration  (outer  iteration  loop)  is  the  Newton-Raphson  iteration,  and  the  bottom-level  iteration 
(inner  iteration  loop)  is  the  PCG  iteration.  The  PCG  iteration  was  applied  to  Eqs.  2. 

At  each  time  step,  three  Newton-Raphson  iterations  and  an  average  of  20  PCG  iterations 
(per  Newton-Raphson  iteration)  were  required  to  obtain  five  significant  digits  of  accuracy  of  the 
response  quantities  obtained  using  Eqs.  2.  The  details  of  the  analysis  are  given  in  Ref.  20. 
Normalized  contour  plots  for  the  displacement  and  velocity  components  at  3.0  msec,  drawn  on 
the  undeformed  middle  surface  of  the  panel,  are  shown  in  Fig.  7.  The  measured  CP  times  and 
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processing  rates  on  a  single  CPU  for  the  different  modules  of  the  proposed  strategy,  are  shown  in 
Fig.  X.  As  can  be  seen  from  Fig.  7,  the  processing  rates  of  most  of  the  modules  was  over  175 
MFLOPS  and  the  processing  rate  for  the  most  time-consuming  module  (module  4  -  incorporation 

I 

of  boundary  conditions  and  matrix  decomposition)  was  274.7  MFLOPS.  Comparison  of  the 
wall-clock  times  obtained  by  the  present  strategy  on  one-,  two-  and  four-CPUs,  with  those  of  the 
direct  analysis  of  the  panel  (with  no  partitioning)  are  given  in  Table  2.  As  can  be  seen  from 
Table  2,  the  use  of  the  present  strategy  on  a  four  CPU  CRAY-YMP  machine  reduces  the  total 
analysis  time  by  over  one  order  of  magnitude,  compared  with  that  required  by  the  direct  analysis 
(on  a  single  processor).  It  is  anticipated  that  the  speedup  ratios  can  be  increased  by  further 
optimizing  the  Fortran  code  used  in  implementing  the  computational  procedure  of  the  partition¬ 
ing  strategy.  The  sustained  processing  rate  would  then  be  increased  to  that  of  the  direct  analysis 
implementation. 


CONCLUDING  REMARKS 

Novel  computational  strategies  are  presented  for  the  analysis  of  large  and  complex  struc¬ 
tures.  The  strategies  are  based  on  generating  the  response  of  the  complex  structure  using  large 
perturbations  from  that  of  a  simpler  structure  (or  a  simpler  mathematical/discrete  model  of  the 
original  structure).  Two  general  approaches  are  described  for  selecting  the  simpler  model  and 
establishing  the  relations  between  the  original  and  simpler  models.  The  two  approaches  are: 
decomposition  or  partitioning  strategy,  and  hierarchical  modeling  strategy.  Two  effective 
partitioning  strategies  are  proposed.  The  first  is  based  on  uncoupling  of  load-carrying 
mechanisms,  and  the  second  is  based  on  symmetry  transformations.  Numerical  results  are 
presented  to  demonstrate  the  effectiveness  of  the  decomposition  strategies. 
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APPENDIX  A  -  PCG  TECHNIQUE 
USED  FOR  THE  SOLUTION  OF  EQS.  3 

For  convenience,  Fx|s.  3  can  be  written  in  the  following  compact  form: 

[l  K  |„  -X|  K  li  |  K  i;,'  |  K  u]  ( 7A  =  I  Q)„  -XIKjx  |  K  |„'  { 0 

where 
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|K|„  = 


Kn  . 

L  •  K22j 


{Q}«,  = 


Qi\ 


.  IK|x 


{zr= 


.  K,2 

LK21  .  . 


The  lour  major  steps  of  the  solution  are: 

Step  1  -  Initialization 

1.  Obtain  an  initial  estimate,  of  (Z)  by  solving  the  equations: 


[KUZ}0  =<Q}0 

Then  obtain  the  vector  {Z}0  using  the  equations: 

I K |0  {Z}„  =  |K|x  {Z}0 

2.  Obtain  the  corresponding  preconditioned  residual,  {y}()  by  solving  the  equations: 

|K|„(y)„  =-|KU((Z}„-{Zj„) 

=  (R(o 

where  {R}„  is  the  initial  residual  vector. 

3.  Set  the  initial  conjugate  search  direction  vectors  {S}„  =  {y}„,  and  obtain  the  vector  -jS}„  by 
solving  the  equations: 

|K|0{S}„  =1K1X{S}„ 

Step  2  -  Line  Search  and  Updating  of  Solution  and  Residual 
For  i=0,  begin  the  iteration  loop  and  do  the  following: 

4.  Compute  the  step  length  at  along  the  search  direction: 

=  (yllWi  _ 

tti  ~  {S}J  {|K|„  { S }  j 

where  superscript  t  denotes  transposition. 
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5.  Update  the  solution: 

{Z}i+I  ={Zh  +  CX;  {  S  }; 

6.  Obtain  the  vector  {Z}j+i  using  the  equations: 

|KUZ}i+,  ={K}x{Z}i+l 

7.  Compute  and  update  the  new  residual  vector  using  the  equations: 

{R}i+i  =  Wi  -Oi|K|0{S}i  +lKJx({Z}i+1-{Z}i) 

Step  3  -  Convergence  Cheek 

S.  If  |R|j+i  <  c  |R|(,  where  IRI  is  the  Euclidean  norm  of  the  residual  vector  and  E  is  a  prescribed 
tolerance,  then  stop,  otherwise  continue. 

Step  4  -  Computation  of  New  Search  Direction  Vector 

9.  Solve  for  the  preconditioned  residual  vector  {y};+i 

|KUy}i+1  ={R}i+1 

10.  Compute  the  orthogonalization  coefficient,  pi+j,  using 

n  _  Wi+iWi+t 
1+1 

1 1 .  Update  the  conjugate  search  direction  vector 

{S}i+1  ={y}i+,  +  Pi+,  {S }j 

12.  Obtain  the  vector  {  S  } ;+1  using  the  equations 

|KUS}i+,  =|K|x{S}i+1 

13.  Step  i  by  1  and  go  to  4. 
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Table  I  -  Number  of  degrees  of  freedom  in  each  partition  for  the  two  decomposition  strategies. 

Cantilevered  composite  shallow  shell  (see  Fig.  1). 


Nonzero 

Degrees 

of 

Freedom 

Uncoupling  of  Load-Carrying 
Mechanisms* 

Symmetry 

Transformations** 

(Z.) 

(Z2) 

(Z,J 

(Z2) 

Stress 

Parameters 

1728 

2880 

2304 

2304 

Generalized 

Displacements 

1200 

1800 

1512 

1488 

*{Z, 1  is  the  vector  of  degrees  of  freedom  associated  with  (N,,  N2,  NI2)  and  (u,,  u2);  and  (Z2|  is 
the  vector  of  degrees  of  freedom  associated  with  (M,,  M2,  MI2,  Q,,  Q2)  and  (w,  <{>,,  <)>.,) 

**{Z, )  and  {ZJ  are  symmetric  and  antisymmetric  components  of  the  response  vector. 


Tabic  2  -  Performance  evaluation  of  the  symmetry  transformation  strategy  on  the  CRAY- 
YMP4/432  at  CRAY  Research,  Inc.,  in  Mendota  Heights,  MN.  UN1COS  5.1.10  operating 
system,  CRT  77  compiler  version  4.0.0. 


Full  Structure 
(optimized  code) 

Partitioned  Structure 
(nearly  optimized  code) 

Number  of 
degrees  of 
freedom 

3  X 1 X  d  i  spl  ace  me  n  t  s 

6144  stresses 

971  displacements 

1536  stresses 

Semi  band  width 
of  equations 

700 

315 

Wall  clock  time 
for  first  ten 
steps,  (sec.) 

171 

58.6  (1  CPU) 

29.7  (2  CPUs) 

16.4  (4  CPUs) 

Sustained  speed 
on  one  CPU, 

M  FLOPS 

277.6 

246 

Speedup  due  to 
multiprocessing 

1.0  (ICPU) 

1.95  (2  CPUs) 

3.45  (4  CPUs) 

Total  speedup 
achieved  by 
partitioning 
strategy 

1.0 

2.92  (ICPU) 

5.76  (2  CPUs) 

10.43  (4  CPUs) 
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Cantilevered  composite  shallow  shell  with  trapezoidal  planform  considered  in  the 
present  study. 

Normalized  contour  plots  for  the  generalized  displacements  evaluated  at  X=() 
(uncoupled  load  carrying  mechanisms)  and  at  A=l.  Cantilevered  composite  shallow 
shell  shown  in  Fig.  1.  Spacing  of  contour  lines  is  0.2  and  dashed  lines  represent 
negative  contours. 

Convergence  of  normal  displacement,  w,  at  point  a  and  total  strain  energy,  U,  ob¬ 
tained  by  using  PCG  technique  in  conjunction  with  Eqs.  2  and  3.  Cantilevered 
composite  shallow  shell  shown  in  Fig.  1.  Iteration  0  corresponds  to  uncoupled 
load-carrying  mechanisms. 

Normalized  contour  plots  for  the  displacement  components  at  A=0  (symmetrized 
structure)  and  at  X=l.  Cantilevered  composite  shallow  shell  shown  in  Fig.  1.  Spacing 
of  contour  lines  is  0.2  and  dashed  lines  represent  negative  contours. 

Convergence  of  normal  displacement,  w,  at  point  a  and  total  strain  energy,  U,  ob¬ 
tained  by  using  PCG  technique  in  conjunction  with  Eqs.  2  and  3.  Cantilevered 
composite  shallow  shell  shown  in  Fig.  1.  Iteration  0  corresponds  to  a  symmetrized 
structure. 

Laminated  anisotropic  composite  panel  with  an  off-center  circular  cutout  used  in  the 
present  study. 

Normalized  contour  plots  for  generalized  displacements  and  velocity  components  at 
t=3.()  Msec.  Laminated  anisotropic  composite  panel  with  an  off-center  circular  cutout 
subjected  to  uniform  normal  loading  po=-50,(MX)  Pa  (see  Fig.  6).  Spacing  of  contour 
lines  is  0.2  and  dashed  lines  represent  negative  contours. 

Measured  CP  times  and  CP  speed  on  one  CPU  of  the  CRAY-YMP4/432.  Laminated 
anisotropic  composite  panel  with  an  off-center  circular  cutout  subjected  to  uniform 
normal  loading  po=-5(),000  Pa  (see  Fig.  6). 
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ABSTRACT 

A  study  is  made  of  a  predictor-corrector  procedure  for  the  accurate  determination  of  the 
temperature  and  heat  flux  distributions  in  thick  multilayered  composite  plates  and  shells.  A 
linear  thnough-the-thickness  temperature  distribution  is  used  in  the  predictor  phase.  The  func¬ 
tional  dependence  of  temperature  on  the  thickness  coordinate  is  then  calculated  a  posteriori  and 
used  in  the  corrector  phase. 

Extensive  numerical  results  are  presented,  for  linear  steady-state  heat  conduction 
problems,  showing  the  effects  of  variation  in  the  geometric  and  lamination  parameters  on  the 
accuracy  of  the  thermal  response  predictions  of  the  predictor-corrector  approach.  Both  antisym¬ 
metrically  laminated  anisotropic  plates,  and  multilayered  orthotropic  cylinders  are  considered. 
The  solutions  are  assumed  to  be  periodic  in  the  surface  coordinates.  For  each  problem  the 
standard  of  comparison  is  taken  to  be  the  analytic  three-dimensional  solution  based  on  treating 
each  layer  as  a  homogeneous  anisotropic  medium.  The  potential  of  the  predictor-corrector 
approach  for  predicting  the  thermal  response  of  multilayered  plates  and  shells  with  complicated 
geometry  is  discussed. 


NOTATION 

c\, ,  c’t  integration  constants  (see  Eqs.  H  and  9) 

h  thickness  of  plate  (or  shell) 

hj_j ,  hj  distance  from  the  middle  surface  to  the  bottom  and  top  surfaces  of  the  jth 

layer,  respectively 

Integrated  thermal  conductivity  coefficients  (see  Eqs.  20) 

„  (i,i=l  to  3;  J=  1,2,3) 

?•<--  - . 

■(t  J>\{  kLL ,  k-rr  thermal  conductivity  coefficients  in  the  direction  of  fibers  and  norma!  to  it. 


1 


ka  (i,i=lto3) 
L, 

L2 

m 

NL 

11 

Q 

Q(,)  (1=1,2) 
cij  (i=l  to  3) 
Rk  (i,i=l  to  3) 


respectively 

thermal  conductivity  coefficients  for  a  typical  layer 
side  length  of  plate  (or  cylinder)  in  X|  direction 
side  length  of  plate  in  X2  direction  (2flr0  for  cylinders) 
longitudinal  wave  number  (in  Xj  direction) 
number  of  layers  of  plate  (or  cylinder) 
circumferential  wave  number 
internal  heat  generation  per  unit  volume 
integrated  heat  generation  (see  Eqs.  17) 
heat  flux  components  in  the  x;  coordinate  directions 
thermal  resistivity  coefficients  for  a  typical  layer 


ro 

SJ.S^  (a,!=l,2) 


T 

T„,Tj 


radius  of  middle  surface  of  cylinder 
integrated  heat  flux  components  (see  Eqs.  17) 
temperature 

temperature  functions  used  in  first-order  heat  conduction  theory  (see  Eq. 


x, ,  x2,  x3 


14) 

orthogonal  coordinate  system 


P 

0 

X  =  1 

=  1  +  x3  /r o 

4.^2,  C 


n, 


tracing  parameter  identifying  anisotropic  terms 

fiber  orientation  angle 
for  plates 
for  cylinders 

dimensionless  coordinates  in  the  \\ ,  x2  and  x3  coordinate  directions  (see 
Fig.  1) 

functional  defined  in  Eq.  5 


Rjt  qjq,  dV  =  thermal  potential  of  structure 


n, 


component  of  n  associated  with  the  in-nlane  heat  flux  components 
q,  and  q2 


2 


n2 


component  of  fl  associated  with  the  transverse  heat  flux  component  q^ 


Ranges  of  indices  (subscripts  and  superscripts) 

I  1,2 

J  1,2,3 

i,  t  1,2,3 
j  I  to  NL 


<x,p  1,2 


INTRODUCTION 

The  expanded  use  of  advanced  composite  materials  in  aerospace,  automotive,  shipbuilding, 
nuclear  and  other  high  tech  industries  has  stimulated  interest  in  the  accurate  prediction  of  the 
thermal  response  of  laminated  anisotropic  plates  and  shells.  A  number  of  approaches  have  been 
developed  for  the  determination  of  the  heat  transfer  characteristics  and  thermal  response  of 
laminated  anisotropic  plates  and  shells.  Some  of  these  approaches  are  similar  to  the  ones  used  in 
predicting  the  mechanical  response  of  composite  plates  and  shells  and  include  microstructural 
models  in  which  the  fibers  and  matrix  are  each  treated  as  homogeneous  isotropic  (or  anisotropic) 
material  (see,  for  example,  Refs.  1,  2  and  3);  and  lamination  theories  in  which  the  multilayered 
composite  is  replaced  by  an  equivalent  homogeneous  anisotropic  material  and  effective  thermal 
conductivities  are  obtained  in  terms  of  the  component  conductivities  (see  Refs.  4,  5  and  6). 
Analytic,  semi-analytic  and  finite  element  solutions  for  the  heat  conduction  three-dimensional 
solids  are  presented  in  Refs.  7,  8,  9,  10,  1 1  and  12.  A  number  of  approaches  have  been  proposed 
for  the  reduction  of  the  three-dimensional  heat  transfer  problem  in  laminated  plates  and  shells  to 
a  two-dimensional  problem.  These  include  global  approximation  models  in  which  the  order  of 
the  governing  differential  equations  is  not  dependent  on  the  number  of  layers  (see,  for  example. 
Ref.  13);  discrete  layer  models  (Refs.  14  and  15)  in  which  the  order  of  the  governing  differential 
equations  depends  on  the  number  of  layers,  integral  transform  techniques  (Ref.  16),  and  operator 
methods  (Ref.  17). 

A  simple  predictor-corrector  procedure  has  been  developed  by  the  authors  in  Refs.  IX,  19 
and  20  for  the  accurate  determination  of  the  global  as  well  as  the  detailed  response  characteristics 
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of  isothermal  plates  and  shells.  The  procedure  is  based  on  using  the  information  obtained  from  a 
simple  two-dimensional  shear  deformation  theory  to  correct  certain  key  elements  of  the  comput¬ 
ational  model  (in  an  inexpensive  postprocessing  mode),  and  hence,  improve  the  response  predic¬ 
tions.  The  present  study  is  an  adaptation  of  the  predictor-corrector  procedure  to  the  heat  transfer 
problem  of  multilayered  composite  plates  and  shells. 

To  sharpen  the  focus  of  the  study,  only  linear  steady-state  heat  conduction  of  composite 
plates  and  cylinders  is  considered.  The  plates  are  assumed  to  be  antisymmetrically  laminated 
with  respect  to  the  middle  plane,  and  the  cylinders  are  constructed  of  orthotropic  layers.  For  both 
the  plates  and  cylinders,  the  thermal  response  quantities  are  assumed  to  be  periodic  in  the  surface 
coordinates. 


2.  MATHEMATICAL  FORMULATION 

Figure  1  shows  the  geometric  characteristics  of  multilayered  plates  and  cylinders  as 
follows:  L|  =  the  length  of  the  plate  (or  cylinder)  in  the  X|  direction;  h  =  total  thickness  of  the 
plate  (or  cylinder);  ru  =  the  radius  of  the  middle  surface  of  the  cylinder;  and  L2  =  the  length  of  the 
plate  in  the  X2  direction  (27cr„  for  the  cylinder).  The  dimensionless  coordinates  £i ,  C  are 
introduced,  where 


(1) 

^  =  rr 2 

(2) 

*3 
=  h 

(3) 

where  A  =  1  for  plates  and  1  +  —  for  cylinders;  and  x?  =  r-r„  for  cylinders. 

*d 

2. 1  Basic  Assumptions 

The  computational  models  are  based  on  the  following  assumptions: 

1.  Each  layer  is  considered  as  a  homogeneous  anisotropic  material  with  the  effective 
thermal  characteristics  obtained  using  micromechanics  equations  (see,  for  example.  Refs.  5  and 
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6). 


2.  Heat  transfer  coefficients  and  thermal  characteristics  of  all  layers  are  independent  of 
temperature. 

3.  Perfect  thermal  contact  exists  between  layers  (i.e.,  the  temperature  and  the  normal  flux 
component  are  continuous  at  layer  interfaces). 

4.  No  heat  flux  is  generated  at  the  internal  surfaces. 

5.  The  thermal  response  is  governed  by  two  sets  of  equations,  namely,  the  heat  flow 
balance  equation  and  the  generalized  Fourier’s  law  of  heat  conduction.  These  equations  are 
given  subsequently  for  multilayered  anisotropic  plates  and  cylinders. 

2.2  Governing  Equations  for  Individual  Layers 

The  governing  three-dimensional  equations  for  each  layer  of  the  plate  (or  shell)  can  be 
written  in  the  following  form  (see,  for  example,  Refs.  21  and  22): 

a)  I  leal  Flow  Balance 

<)|  qi  +  ^2 M2  +  \  Ms)  ~  Q  =  0  (4) 

and 

b)  Generalized  Fourier’s  Law 


^11  P  *12  P  k\i 

•3|T| 

M2 

A:  22  P^23 

< 

*T 

Uu. 

_symm  £33  _ 

UjTj 

or,  the  inverse  relations 

r  R„  pR12  0R|.r 

1*22  P  1*2.1 

.symm  R33  J 


re _ r 

v  2  1 

J}' T> 

V 

1 

1 

(6) 


where  q;  (i=l  to  3)  are  the  flux  components  in  the  X|  coordinate  directions;  T  is  the  temperature; 
k\l=kl\  (i,t=l  to  3)  are  the  thermal  conductivity  coefficients  of  the  material  of  the  layer;  R„  are  the 
thermal  resistivity  coefficients;  Q  is  the  internal  heal  generation  per  unit  volume;  d,  h  d/dx,;  and 

P=()  and  1  for  orthotropic  and  anisotropic  layers,  respectively. 


Equations  4  and  6  arc  the  stationary  conditions  of  the  following  two-field  functional: 


n,  (Mi  ,T)  =  J 


ili  d,  T  +  q,  q,  Ra  +  QT 


dxj  dx2  dxy 


(7) 


where  d j  =  (i  =  1  to  3),  and  a  repeated  index  i  or  t  in  the  same  term  denotes  summation  over 


its  full  range  (1  to  3). 

2.3  Reduction  to  a  Two-Dimensional  Problem 

In  the  present  study,  the  reduction  of  the  three-dimensional  heat  conduction  problem 
described  in  Eqs.  4  and  5  (or  6)  to  a  two-dimensional  problem  is  accomplished  by  assuming  a 
linear  temperature  variation  in  the  thickness  coordinate;  and  replacing  the  laminated  plate  (or 
shell)  by  an  equivalent  single-layer  plate  (or  shell).  The  governing  equations  of  the  resulting 
first-order  heat  conduction  theory  are  obtained  by  integrating  Eqs.  4  and  5  in  the  thickness 
direction.  Because  of  the  discontinuity  of  the  thermal  conductivity  coefficients,  the  integration  is 
performed  layer-by-layer.  The  fundamental  equations  of  the  first-order  heat  conduction  theory 
are  given  in  the  Appendix. 


3.  PREDICTOR-CORRECTOR  PROCEDURE 

3. 1  Basic  Idea  of  the  Procedure 

The  predictor-corrector  procedure  used  in  the  present  study  is  an  iterational  process  in 
which  the  thermal  response  obtained  in  the  first  (predictor)  phase  of  the  analysis  is  used  to 
correct  the  temperature  distribution  and,  hence,  improve  the  response  predictions.  Numerical 
experiments  have  shown  that  only  one  or  two  iterations  are  needed  (in  the  corrector  phase)  to 
obtain  highly  accurate  response  predictions.  The  application  of  the  procedure  to  multilayered 
composite  plates  and  shells  is  described  subsequently.  The  superscript  o  refers  to  the  predictions 
of  the  first-order  heat  conduction  theory;  and  a  bar  (-)  over  a  symbol  refers  to  the  thermal  re¬ 
sponse  quantities  obtained  by  using  three-dimensional  equatiops. 

3.2  Predictor  Phase 

The  first-order  theory  outlined  in  the  previous  section  is  used  to  evaluate  through-the- 
thickness  temperature,  T’,  and  flux  components  q"  and  q\’.  Then  the  flux  component  <|(  by 
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integrating  Eqs.  1  in  the  thickness  direction  as  follows: 

*3 

q;>  =  -v-  J  c]Y  +  (h q,  -O^dx^  +  c,j  (X) 

-li/2 

The  temperature  distribution  in  the  thickness  direction,  T,  is  obtained  by  integrating  the 
third  equation  of  Eqs.  3  as  follows: 

\3 

T=  J  ["p^Rnq'i’  +R23M2)  +  R33  43]dx3  +  ct  (tJ) 

-i  1/2 

In  Eqs.  X  and  9,  ct|  and  cj  are  integration  constants  obtained  from  the  conditions  at  the  outer 
surfaces  of  the  laminate. 

Note  that  because  of  the  discontinuity  of  c]|  and  q2  at  layer  interfaces,  the  integrations  in 
Eqs.  8  and  9  are  performed  in  a  piecewise  manner  (layer  by  layer). 

3.3  Corrector  Phase 

The  calculation  of  the  corrected  thermal  response  of  the  plate  (or  shell)  may  be  con¬ 
veniently  divided  into  three  steps:  namely,  1)  generation  of  coordinate  (basis)  temperature 
functions;  2)  computation  of  amplitudes  of  the  coordinate  functions,  and  3)  evaluation  of  the 
corrected  through-the-thickness  flux  components  using  three-dimensional  continuum  equations. 
The  first  two  steps  are  described  subsequently. 

1.  Generation  of  Coordinate  (Basis)  Temperature  Functions.  The  temperature  T  is 
decomposed  into  symmetric  and  antisymmetric  functions  of  the  thickness  coordinate  x>  Each  of 
the  symmetric  and  antisymmetric  functions  is  further  subdivided  into  an  initial  and  correction 
function.  The  initial  functions  are  associated  with  the  temperature  distribution  T°,  obtained  in  the 
predictor  phase.  The  correction  functions  are  nonlinear  functions  in  x^  and  are  associated  witli 
the  difference  T-T\  Since  the  initial  functions  satisfy  the  temperature  boundary  conditions  at 

l' 

the  outer  surfaces  of  the  plate  (or  shell),  the  correction  functions  vanish  on  these  surfaces. 

2.  Computation  of  the  Amplitudes  of  the  Coordinate  Functions.  The  resulting  four 
symmetric/antisymmetric  functions  associated  with  the  temperature  arc  now  chosen  as  eoor- 
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dinate  (or  basis  functions),  and  the  temperature  is  expressed  as  a  linear  combination  of  the  four 
functions,  with  unknown  parameters  (representing  the  amplitudes  of  the  coordinate  functions  (or 
temperature  modes).  The  four  unknown  parameters  are  obtained  by  using  a  direct  variational 
procedure. 

Note  that  further  improvement  can  be  obtained  by  performing  additional  correction  steps. 
This  is  accomplished  by  replacing  T’.q"  and  in  Eqs.  X  and  9  by  the  response  quantities 

obtained  in  the  preceding  corrector  step  and  repeating  the  corrector  phase. 

4.  NUMERICAL  STUDIES 

To  assess  the  accuracy  and  effectiveness  of  the  predictor-corrector  computational  proce¬ 
dure,  a  large  number  of  heat-conduction  problems  of  multilayered  composite  plates  and  cylinders 
have  been  solved  by  these  techniques.  The  composite  plates  considered  in  the  present  study  are 
square  laminates  with  L!  =L2  =  1.0,  and  have  antisymmetric  lamination  with  respect  to  the  middle 
plane.  The  composite  cylinders  considered  are  closed  circular  cylinders  with  orthotropic  layers, 
and  with  L|  =r„  =  1.0.  The  fibers  of  the  different  layers  alternate  between  the  circumferential  and 

longitudinal  directions,  with  the  fibers  of  the  top  layers  running  in  the  circumferential  direction. 
The  temperature  boundary  conditions  at  the  top  and  bottom  surfaces,  as  well  as  the  thermal 
response,  are  periodic  in  X|  and  X2  with  periods  2L\  and  2Lj  for  plates  (2L|  and  L2  for 
cylinders). 

For  both  the  plates  and  the  cylinders,  the  X3  coordinate  is  a  principal  material  direction. 
Moreover,  for  the  antisymmetrically  laminated  plates,  the  thermal  conductivity  coefficients  for  a 
pair  of  layers  j+  and  j‘  which  are  symmetrically  situated  with  respect  to  the  middle  plane  satisfy 
the  following  relations: 


II 

(i=l  to  3  and  i 

is  not  summed) 

£  c: 

— 

II 

£  O 
vj 

II 

(a=l,2) 

(10) 

;.U0  _  k«  ) 

*  12  ~  *  12 

1 

with  similar  relations  for  the  thermal  resistivity  coefficients. 

Equations  10  are  satisfied  for  angle-ply  antisymmetric  laminates,  cross-ply  symmetric 
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laminates,  and  combinations  of  the  two.  In  the  case  of  symmetric  or  unsymmetric  cross-ply 
laminates,  k  12  =0  and  the  matrix  of  thermal  conductivity  coefficients  is  diagonal.  The  same  is 
true  for  the  matrix  of  thermal  resistivity  coefficients. 

The  material  characteristics  of  the  individual  layers  are  taken  to  be  those  typical  of  high- 
modulus  fibrous  composites,  namely: 

*ll=5-0  ,  &tt=0-5 

where  subscript  L  refers  to  the  direction  of  fibers  and  subscript  T  refers  to  the  transverse  direc¬ 
tion.  For  plates  the  symmetric  and  antisymmetric  (in  the  thickness  direction)  components  of  the 
temperature  were  assumed  to  be  of  the  form: 

Ts  sin  rt^i  sin  7r^2  and  (11) 

Ta  cos  cos  jc£,2  (12) 


where  subscripts  s  and  a  refer  to  the  symmetric  and  antisymmetric  components  (in  the  thickness 
direction).  The  values  of  Ts  and  Tu  at  the  top  surface  of  the  plate  are  given  by 
Ts  =  250  and  T„  =  50.  Note  that  the  trigonometric  functions  used  in  describing  the  surface 

distributions  of  the  symmetric  and  antisymmetric  components  of  the  temperature  are  different. 

For  cylinders,  the  temperature  distribution  was  assumed  to  be  of  the  form: 


) 


cos  2nn  £,2 


(13) 


where  Tn=  300  and  T;  =  200.  Subscripts  o  and  i  refer  to  the  outer  (top)  and  inner  (bottom) 


surfaces  of  the  cylinder,  respectively.  For  each  problem,  the  solutions  obtained  by  the  four 
models  described  in  Table  1  were  compared  with  the  analytic  three-dimensional  continuum 
solutions. 
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Table  1  -  Characteristics  of  the  models  used  in  the  numerical  studies 


Model 

Predictor  Phase 

Number  of  iterations  in 
the  corrector 

1 

First-order  heat  conduction  theory 

0 

2 

(based  on  assuming  a  linear  temperature 

1 

2A 

variation  in  the  thickness  coordinate) 

2 

2B 

3 

For  plates,  three  parameters  were  varied,  namely:  the  thickness  ratio  of  the  plate,  h/L| ;  the 
number  of  layers,  NL;  and  the  fiber  orientation  angle  of  the  individual  layers,  0.  The  thickness 
ratio  was  varied  between  0.01  and  0.5;  the  number  of  layers  was  varied  between  2  and  20;  and  0 
was  varied  between  0°  and  45°.  For  cylinders,  three  parameters  were  varied,  namely:  the 
thickness  ratio  h/r();  the  number  of  layers,  NL;  and  the  circumferential  wave  number,  n.  The 
longitudinal  wave  number  was  selected  to  be  1,  and  the  length-to-radius  ratio,  L/r,„  was  selected 
to  be  1.0.  The  number  of  layers  was  varied  between  2  and  20;  h/r0  between  0.0 1  and  0.5,  and  n 
between  0  and  10. 

As  a  step  towards  establishing  the  range  of  validity  of  the  predictor-corrector  procedure 
and  the  number  of  iterations  required  in  the  corrector  phase,  the  thermal  potential  of  the  structure, 

n  =  J  Rii  Mi  Mi  dV  (i,i  =  1  to  3);  was  decomposed  into  two  components:  n i  associated  with 

qt  and  c)2;  and  Il2  associated  with  q3  (112  =  \  2  R33  M3  M3  dV).  The  total  thermal  potential  of 

the  structure  fl  =  FI  j  +  The  assessment  of  the  predictor-corrector  procedure  included  both 
the  global  thermal  response  characteristics,  n,  n j  and  112,  as  weH  as  detailed  temperature  and 
flux  distributions  in  the  thickness  direction. 

Typical  results  are  given  in  Figs.  2,  3  and  4  for  the  antisymmetrically  laminated  plates,  and 
in  Figs.  5,  6  and  7  for  the  multilayered  orthotropic  cylinders.  The  effects  of  variation  of  the  two 
parameters,  h/L| ,  NL  for  plates  (and  h/L,  and  n  for  cylinders)  on  the  thermal  potential  compo- 
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nents  I-!]  and  riot  obtained  by  the  three-dimensional  continuum  model  are  depicted  in  Fig.  2 
for  plates  and  in  Fig.  5  for  cylinders. 

An  indication  of  the  accuracy  of  the  thermal  potential  components  n |  and  n2,  obtained  by 
the  predictor-corrector  procedure,  is  given  in  Figs.  3  and  6.  Figures  4  and  7  give  an  indication  of 
the  accuracy  of  the  temperature  and  heat  flux  distributions  in  the  thickness  direction.  In  Fig.  4 
both  the  symmetric  and  antisymmetric  parts  of  the  thermal  response  quantities  (with  respect  to 
the  middle  plane)  are  shown.  Note  that  since  the  symmetric  and  antisymmetric  components  of 
each  thermal  response  quantity  are  multiplied  by  different  trigonometric  functions  in  X|  and  x2, 

the  value  of  the  response  quantity  is  a  linear  combination  of  the  two  components. 

An  examination  of  Figs.  2  to  7  reveals: 

1.  The  ratio  of  the  thermal  potential  component  I~l]  /fl  increases  rapidly  with  the  increase 
in  h/L|  from  0.01  up  to  h/Li=0.1  for  plates  (and  h/Lj  from  0.01  to  0.15  for  cylinders  with  n  <  2). 
The  decrease  in  n |  /n  is  associated  with  a  sharp  decrease  in  n2/n.  Further  increases  in  h/L|  are 
associated  with  a  slow  decrease  in  n i  /n  and  a  slow  increase  in  n2/n  (see  Figs.  2  and  5). 

2.  The  ratios  n!/n  and  n2/n  are  somewhat  insensitive  to  variations  in  the  number  of 
layers  NL.  For  plates  these  ratios  are  also  somewhat  insensitive  to  the  fiber  orientation  angle  0 
(results  not  shown). 

3.  For  antisymmetrically  laminated  plates,  the  accuracy  of  the  predictions  of  the  first-order 
heat-conduction  theory.  Model  1 ,  deteriorates  rapidly  as  h/L|  increases.  This  is  true  for  both  the 

global  as  well  as  the  detailed  thermal  response  characteristics.  For  plates  with  NL=1(),  0=45”, 
h/L|=().  1,  the  error  in  n  predicted  by  Model  1  was  7.4%  and  increased  to  165%  for  h/L|=0.3  (see 
Fig.  3). 

4.  The  accuracy  of  the  predictions  of  the  first-order  theory  and  its  range  of  validity  can  be 
significantly  improved  by  using  one  or  two  iterations  in  the  corrector  phase.  Models  2  and  2A. 
even  for  very  thick  plates  and  shells.  As  an  example  to  this,  for  plates  with  h/L|=0.5,  the  error  in 
11  was  249%.  The  corresponding  errors  after  one  and  two  iterations  in  the  corrector  phase  were 
only  3.2%  and  0.08%,  respectively  (see  Fig.  3). 

5.  For  orthotropic  cylinders  the  accuracy  of  the  predictions  of  Model  1  deterioriates 
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rapidly  as  the  circumferential  wave  number,  n,  increases.  As  an  example  to  this,  for  cylinders 
with  h/r  =0.5,  NL=1(),  n=2,  the  error  in  n  predicted  by  Model  1  w-as  220%,  and  increased  to 

378%  for  n=6.  The  use  of  one  and  two  iterations  in  the  corrector  phase  (Models  2  and  2A) 

I 

reduced  the  errors  for  n=6  to  9.8%  and  0.64%,  respectively  (see  Fig.  6). 

6.  The  thickness  distributions  of  temperature  and  flux  components  obtained  by  the 
predictor-corrector  approach  are  highly  accurate.  For  plates  and  cylinders  with  h/L|  <  0.5  and  n 
<  5,  (lie  distributions  of  the  response  quantities  after  two  iterations  (Model  2A)  in  the  corrector 
phase  were  almost  indistinguishable  from  the  exact  three-dimensional  solutions  (see,  for  ex¬ 
ample,  Figs.  4  and  7). 

7.  Numerical  experiments  (not  presented  herein)  have  shown  that  the  use  of  higher-order 
heat  conduction  theories  for  multilayered  composites  based  on  a  global  cubic  variation,  or 
higher-degree  polynomial  approximation,  for  the  temperature  through  the  thickness  results  in 
highly  accurate  temperature,  as  well  as  heat  flux  components  q i  and  i]2-  However,  the  transverse 
heat  flux  component  q3  obtained  by  these  theories  is  not  as  accurate. 

5.  POTENTIAL  OF  THE  PREDICTOR-CORRECTOR  PROCEDURE 

The  predictor-corrector  procedure  appears  to  have  high  potential  for  the  accurate  predic¬ 
tion  of  the  thermal  response  of  multilayered  composite  plates  and  shells.  The  numerical  studies 
conducted  for  antisymmetrically  laminated  anisotropic  plates  and  simply  supported  orlhotropic 
cylinders  demonstrated  the  accuracy  and  effectiveness  of  the  predictor-corrector  procedure.  In 
particular,  the  following  two  points  are  worth  mentioning: 

1.  The  predictor-corrector  procedure  can  be  applied,  in  conjunction  with  finite  element 
models,  to  the  analysis  of  anisotropic  plates  and  shells  with  arbitrary  geometry.  The  calculation 
of  the  heat  flux  component  in  the  transverse  direction,  and  the  correction  phase  (including  the 
calculation  of  the  transverse  temperature  distribution)  can  be  performed  on  the  element  level  for 
selected  elements  (in  the  critical  regions  of  the  plate  and  shell, models). 

2.  Although  any  of  the  two-dimensional  heat  conduction  models  can  be  used  in  the  first 
(predictor)  phase  of  the  predictor-corrector  procedure,  the  first-order  theory  involves  fewer 
temperature  parameters  than  higher-order  theories. 
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6.  CONCLUDING  REMARKS 


A  study  is  made  of  a  predictor-corrector  procedure  for  the  accurate  determination  of  the 
temperature  and  heat  flux  distributions  in  thick  multilayered  composite  plates  and  shells.  A 
first-order  heat  conduction  theory,  based  on  linear  through-the-thickness  temperature  distribu¬ 
tion,  is  used  in  the  predictor  phase.  The  functional  dependence  of  temperature  on  the  thickness 
coordinate  is  then  calculated  a  posteriori  and  used  in  the  corrector  phase. 

Extensive  numerical  results  are  presented,  for  linear  steady-state  heat  conduction 
problems,  showing  the  effects  of  variation  in  the  geometric  and  lamination  parameters  on  the 
accuracy  of  the  thermal  response  quantities  obtained  by  the  predictor-corrector  approach.  Both 
antisymmetrically  laminated  anisotropic  plates,  and  multilayered  orthotropic  cylinders  are 
considered.  The  solutions  are  assumed  to  be  periodic  in  the  surface  coordinates,  and  for  each 
problem  the  standard  of  comparison  is  taken  to  be  the  analytic  three-dimensional  solution  based 
on  treating  each  layer  as  a  homogeneous  anisotropic  medium. 

The  numerical  results  clearly  demonstrate  the  effectiveness  of  the  predictor-corrector 
procedures  for  the  accurate  determination  of  the  global  as  well  as  detailed  thermal  response 
characteristics  of  multilayered  plates  and  shells.  The  accuracy  of  the  response  quantities  ob¬ 
tained  in  the  first  (predictor)  phase  for  thick  laminates  (h/L|  >  0.1)  may  be  unacceptable. 
However,  the  use  of  one  or  two  iterations  in  the  corrector  phase  improves  the  predictions  sub¬ 
stantially  and  results  in  highly  accurate  distributions  of  temperature  and  heat  flux  components 
through  the  thickness. 
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APPENDIX  -  FUNDAMENTAL  EQUATIONS  OF  THE  FIRST-ORDER 
HEAT  CONDUCTION  THEORY  USED  IN  THE  PRESENT  STUDY 

The  fundamental  equations  of  the  two-dimensional  first-order  heat  conduction  theory  used 

in  the  present  study  are  given  in  this  appendix. 

Temperature  Distribution 

The  temperature  is  assumed  to  have  a  linear  variation  in  the  thickness  coordinate  as  follows: 
T  =  T„  +x3T,  (14) 

where  T„  and  T|  are  functions  of  X|  and  x2. 

Heat  Balance  Equations 

The  two  global  heat  balance  equations  can  be  written  in  terms  of  the  integrated  heat  flux 
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components  as  follows: 

da  S(J>  -  Q(l)  +  (\q3)|  _  ^  -  (kq3)  |  _  =0 

*  'aj— 

s;;>  -  s< 1 1 5  -  q(2)  +  ( x3  Xq3)  i  -  ( x3  Xq3)  i  =  o 

'x3  =h/2  'x3  =-li/2 

where 


f  s">' 

NL  hi 

rv  _  vd-i)7 

Y(<x)  C1(X  x3 

si" 

r=$  j  \ 

Xq3 

.  0">, 

J  hH 

s.  *-Qx3-I)  J 

(15) 

(16) 


(17) 


with  Y(i)  -X  and  7(2)  =  1,  a  =  1,2  and  1  =  1,2.  Note  that  a  repeated  index  a  in  Eqs.  15  and  16 
denotes  summation  over  the  range  1,2;  and  a  is  not  summed  in  Eqs.  17. 

Generalized  Fourier’s  Law 


sS,=-K«  OpT.-KST,  -K^apT, 
S',”  ^-Kji’apTo-K'l’Ti  -Kfj’apT, 


(IS) 

(19) 


where 

NL  hi 

Kj?  =  J  v(i)  V(1)  *ii  dx3  (20) 

J=1  V. 

with  V(|)  =  V(3)  =  VaT  ,  V(2)  =  1/VX  ;  J  =  1,2,3;  i,i  =  1,2,3;  and  a,P  =  1,2.  Note  that  a  repeated 

index  p  in  Eqs.  18  and  19  denotes  summation  over  the  range  1,2;  and  i,i  are  not  summed  in  Eqs. 

20. 

Governing  Equations 

The  two  governing  differential  equations  in  the  temperature  functions  T„  and  T(  are  oh- 

I* 

tained  by  replacing  the  quantities  in  Eqs.  14  and  15  by  their  expressions  in  terms  of 

T„  and  T|  (Eqs.  19  and  20). 

Trigonometric  Functions  for  the  Temperature  and  Integrated  Heal  Flux  Components 
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The  temperature  and  integrated  heat  flux  components  are  assumed  in  the  present  study  to  be 
periodic  in  both  the  \|  and  x2  directions  with  periods  of  2L|  and  2L2  for  plates  (and  2Lj ,  L2  for 
cylinders).  The  following  products  of  trigonometric  functions  are  used  for  the  different  response 

i 

quantities: 

Antisymmetrically  laminated  plates: 


T„  =T„  sin  mit^i  sin  n7t£2 

mn  J 

(21) 

(T"s-',)=(T|™’s’l)c'osm^ 

cos  ntt^2 

(22) 

(s<’\  s<2>  Ws^  ,s(22)  ) 

V  /  \  mn  ^mn/ 

cos  twr^i 

cos  n^£,2 

(23) 

=  ("s(21}  ,s(2)  ) 

V  *  /  \  Lmn  1  mn  / 

sin  mnS,] 

cos  nn^2 

(24) 

Multilayered  orthotropic  cylinders: 

(T,..TI.S<,»)=(T„mn,T|m„.S«'»in) 

cos  mit^i  cos  2jtn^2 

(25) 

sin  nm^i 

cos  2jtn^2 

(2b) 

OM)  -(S^.S®) 

cos  mJtE,[ 

sin  27tn^2 

(27) 

The  governing  equations  can  be  reduced  to  a  set  of  algebraic  equations  through  the  use  of  the 
trigonometric  functions,  Eqs.  21  to  27.  Note  that  the  trigonometric  functions,  Eqs.  21  to  27, 
provide  exact  solution  for  the  governing  differential  equations  of  antisymmetrically  laminated 
plates  and  multilayered  orthotropic  cylinders. 
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Figure  I  -  Laminated  composite  plates  and  cylinders  used  in  the  present  study. 

Figure  2  -  Effect  of  thickness  ratio,  h/Lj,  and  number  of  layers,  NL,  on  components  of  the 
thermal  potential  obtained  by  three-dimensional  continuum  model.  Antisymmetri¬ 
cally  laminated  plates  with  prescribed  surface  temperatures  Ts  sin  7it,|  sin  nt,2  and 

T„  cos  rcS,!  cos  nt,2.  At  the  top  surface  Ts  =  250,  T„=  50,  9=45°. 

Figure  3  -  Effect  of  thickness  ratio,  h/Lj,  on  the  accuracy  of  the  total  thermal  potential,  fl,  and 
the  thermal  potential  components,  I~Ii ,  112,  obtained  by  the  predictor-corrector  proce¬ 
dure.  Antisymmetrically  laminated  plates  with  prescribed  surface  temperatures 
Ts  sin  sin  nt,2  and  Ta  cos  cos  nt,2-  At  the  top  surface 

Ts=  250,  Tu  =  50,  0=45°,  and  NL=10. 

Figure  4  -  Accuracy  of  temperature  and  heat  flux  components  obtained  by  predictor-corrector 
procedures.  Antisymmetrically  laminated  plates  with  prescribed  surface  temperatures 

Ts  sin  7t£,i  sin  7t^2  and  Tu  cos  i  cos  7t^2.  At  the  top  surface 

T,  =  250,  Ta=  50,  9=45°,  h/L,=0.5,  and  NL=10. 

Figure  5  -  Effect  of  thickness  ratio,  h/Lj ,  and  circumferential  wave  number,  n,  on  thermal 
potential  components,  obtained  by  the  three-dimensional  continuum  model.  Multi¬ 
layered  composite  cylinders  with  prescribed  surface  temperatures 
T„  =  T„  cos  jc£i  cos  27tnt,2»  T;  =T;  cos  cos  2nn^2>  T,>=  300,  Tj  =  200,  L/r„ 

=  1.0,  NL=10. 

Figure  6  -  Effect  of  circumferential  wave  number,  n,  on  the  accuracy  of  the  total  potential,  11. 

and  thermal  potential  components,  11)  ,  n2,  obtained  by  the  predictor-corrector 
procedures.  Multilayered  composite  cylinders  with  prescribed  surface  temperatures 
T„  =  T„  cos  7t£, |  cos2nn^2,  T;=TjCosn^|  cos  2nn£,->;  T„  =  300, T;  =  200,  L/r„ 

=  1.0,  NL=  10,  h/r  =0.5. 


Figure  7  -  Accuracy  of  temperature  and  heat  flux  components  obtained  by  predictor-corrector 
procedure.  Multilayered  composite  cylinders  with  prescribed  surface  temperatures 
T0  =  T„  cos  7r£,i  cos  27m£,2,  T,  =T;  cos  cos  Inn^  T„=  300,  T,=  200,  L/r„ 

=  1.0,  h/r  =0.5,  NL=1(),  and  n=l. 


